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The hadron string dynamics (HSD) model is generalized to include the creation and evolution 
of retarded electromagnetic fields as well as the influence of the magnetic and electric fields on 
the quasiparticle propagation. The time-space structure of the fields is analyzed in detail for non- 
central Au-|-Au collisions at ^/sNN =200 GeV. It is shown that the created magnetic field is highly 
inhomogeneous but in the central region of the overlapping nuclei it changes relatively weakly 
in the transverse direction. For the impact parameter b =10 fm the maximal magnetic field - 
perpendicularly to the reaction plane - is obtained of order eBy/m^ ~5 for a very short time ~ 
0.2 fm/c, which roughly corresponds to the time of a maximal overlap of the colliding nuclei. We 
find that at any time the location of the mgiximum in the eBy distribution correlates with that of 
the energy density of the created particles. In contrast, the electric field distribution, being also 
highly inhomogeneous, has a minimum in the center of the overlap region. Furthermore, the field 
characteristics are presented as a function of the collision energy and the centrality of the collisions. 
To explore the effect of the back-reaction of the fields on hadronic observables a comparison of HSD 
results with and without fields is exemplified. Our actual calculations show no noticeable influence 
of the electromagnetic fields - created in heavy-ion collisions - on the effect of the electric charge 
separation with respect to the reaction plane. 

PACS numbers: 25.75.-q, 25.75.Ag 



I. INTRODUCTION 



The study of the properties of QCD matter in the pres- 
ence of strong uniform magnetic fields has attracted much 
attention during recent years due to several remarkable 
observations. They include such a universal phenomenon 
as the magnetic catalysis [J-ll] in which the magnetic field 
acts as a strong catalyst of dynamical flavor symmetry 
breaking which might lead to the generation of a dy- 
namical fermion mass. Furthermore, in dense QCD mat- 
ter in the presence of an external magnetic field and/or 
topoloigical defects, a spontaneous creation of axial cur- 
rents |j| may happen. Moreover, at finite baryon density, 
due to a response of the QCD ground state to a strong 
magnetic field, a metastable object, the 7r° domain wall 



(or "Goldstone current state" in quark matter), could 
appear which energetically may be more favorable than 
nuclear matter at the same density The presence of 
a magnetic field also favors the formation of spatially in- 
homogeneous spiral-like quark condensate configurations 
at low temperatures and non-zero chemical potentials 0] . 
The influence of a constant magnetic field on possible 
color-superconducting phases (the color Meissner effect) 
has also actively been discussed Q- A clarification of 
such phenomena by experimental observations requires 
e.g. the production of QCD matter in relativistic heavy- 
ion collisions where in non-central reactions strong elec- 
tromagnetic fields are created by the charged four-current 
of the spectators. 

In this general respect exact solutions of the quantum 
field equations of motion are of special interest. The lat- 



2 



tcr supply us with microscopic insight for problems of rel- 
ativistic charged particle motion in electromagnetic fields 
of terrestrial experimental devices as well as in astro- 
physics and cosmology. In particular, they apply to the 
development of the quantum theory of synchrotron radi- 
ation |8| and also for the description of interacting parti- 
cles, including electrons and neutrinos, especially in mat- 
ter in the presence of external electromagnetic fields . 

More specifically, the quark-hadron and chiral symme- 
try restoration transitions might be dramatically modi- 
fied in the presence of a strong magnetic field. We recall 
that without a magnetic field a crossover at T = Tc wlGO 
MeV is realized at vanishing baryon chemical potential 
liq. The presence of a strong magnetic background turns 
this picture within a linear sigma model into a weak first 
order phase transition [l^, [HI whereas in the Nasmbu- 
Jona-Lasinio N JL model the crossover remains [l^, [13 ■ 
This contradiction is reconciled by noting a crucial role 
of of the vacuum contribution from quarks. The vacuum 
contribution seems to soften the order of phase transition: 
the first order phase transition - which would be realized 
in the absence of the vacuum contribution - becomes a 
smooth crossover if the system with vacuum quark loops 
included 

Within the Polyakov-Nambu-Jona-Lasinio (PNJL 
model) it was shown that the external field acts as a cat- 
alyzer for dynamical symmetry breaking which increases 
the critical temperature with increasing st reng th of the 
magnetic field and sharpens the transition |15| in agree- 
ment with lattice QCD calculations [l^. Very recently, 
an astonishing feature has been demonstrated in effective 
models in that lines of the finite temperature deconfine- 
mcnt and chiral transitions can split in a strong magnetic 
background [13 ■ We also mention that the presence 
of a strong constant magnetic field modified also the na- 
ture of the electroweak phase transition in the evolution 
of the universe at its early stages (Tsj . 

In the last few years, particular attention has been paid 
to the chiral magnetic effect (CME) closely related to a 
possible local V and CV symmetry violation in strong in- 
teractions as su gge sted in Ref [19j and widely discussed 
in Refs. [l^ [20l423| . This effect originates from the ex- 
istence of nontrivial topological configurations of gauge 
fields and their interplay with the chiral anomaly which 
results in an asymmetry between left- and right-handed 
quarks. Such a chiral asymmetry, when coupled to a 
strong magnetic field as created by colliding nuclei per- 
pendicularly to the reaction plane, induces an electric 
charge current along the direction of a magnetic field 
thereby separating particles of opposite charges with re- 
spect to the reaction plane. Thus, such topological ef- 
fects in QCD might be observed in heavy-ion collisions 
directly in the presence of very intense external electro- 
magnetic fields due to the "Chiral Magnetic Effect" as 
a manifestation of the spontaneous violation of the CV 
symmetry. Recently, topological charge fluctuations and 
possible CME have been confirmed by QCD lattice cal- 
culations in quenched SU{2) gauge theory [^J] and in 



QCD-I-QED with dynamical 2 -|- 1 quark flavors (25j . 

One should note that in contrast to all the cases men- 
tioned above, the magnetic fleld in the CME is not con- 
stant and acts only during a very short time. The maxi- 
mal strength of the electromagnetic field eBy created in 
relativistic heavy-ion collisions is very high eBy ~5 rn^ 
but its duration is t ^0.2 fm/c for Au-hAu at ^Jsnn =200 
GeV collisions and impact parameter h =10 fm, as was 
estimated dynamically in [2l|, [2^, [23| ■ 

It is remarkable that the proposed CME observable, 
i.e. the electric charge asymmetry of produced parti- 
cles with respect to the reaction plane, has been recently 
measured by the STAR Collaboration [Ullil. Qualita- 
tively the experimental results agree with the magnitude 
and gross features of the theoretical predictions for local 
7^— odd violation in heavy-ion collisions. The observed 
signal can not be described by the background models 
used in j28l. [29j. however, alternative mechanisms result- 
ing in a similar charge separation effect are not fully ex- 
cluded (see e.g. Refs. j3(]| - l35| ). This issue is under in- 
tensive debate now. 

The aim of this paper is to study the space-time evo- 
lution of (electro-) magnetic fields formed in relativis- 
tic heavy-ion collisions. The Hadron String Dynamics 
(HSD) transport code [s^ is used as a basis of our con- 
siderations. In contrast with the UrQMD model with- 
out including electromagnetic fields (previously used for 
the analysis jl^, H^), the HSD model corresponds to 
Kadanoff-Baym rather than Boltzman kinetic equations 
and treats the nuclear collisions in terms of quasiparticles 
with a finite width. In our approach the dynamical for- 
mation of the electromagnetic field, its evolution during 
a collision and influence on the quasiparticle dynamics as 
well as the interplay of the created magnetic and electric 
fields and back-reaction effects are included simultane- 
ously. 

The article is organized as follows: In Sect. II the 
model is presented and it is shown how the formation 
of the electromagnetic field and particle propagation in 
this background can be implemented in the HSD trans- 
port code. In Sect. HI the space-time structure of the 
formed electromagnetic fields, their correlations with the 
energy density of produced particles, the correlation be- 
tween the magnetic and electric fields and some estimates 
for the dependence of these characteristics on collision en- 
ergy and impact parameter are investigated. A compari- 
son of some observables, which are calculated within the 
HSD model with and without the electromagnetic fields, 
is presented in Sect. IV. The results are summarized in 
Sect. V. 



II. ELECTROMAGNETIC FIELD EVOLUTION 
WITHIN THE HSD MODEL 

To describe a collision of heavy ions let us start from 
the relativistic Boltzmann equation for an on-shell phase- 



3 



space distribution function / = f(x,p) 

where C[f] is the collision integral and x,p are the 4- 
coordinate and 4-momentum of a particle. A background 
electromagnetic field may be taken into account by in- 
cluding an electromagnetic tensor into Eq.dT]) as 



d 



(2) 



where 



sFni, — dn,Au — dp Art 



with the electromagnetic 4-vector potential = 
{$, A}. Note that the left-hand side of Eq.© IS gauge 
invariant. 

Let US rewrite it terms of components of and 
generalize it to the case of a set of particles moving in a 
nuclear potential U. Then the equation for a test particle 
characterized by the distribution function / = /(r,p,i) 
may be presented as follows: 

OA 



Vr C/ + eV<I> + 

Ccouif, fl, ■■Jn) 



ev X (V X A) ) Vp } f 

(4) 



B = V X A 



E 



(5) 



The strength of the magnetic B and electric E fields is, 
respectively, 

OA 
'dt ■ 

One should note that the electromagnetic field generated 
by moving nuclei may be considered as an external field: 
the value of the electromagnetic field at a given point 
is determined by the global charge distribution of collid- 
ing nuclei and thus, in good approximation, independent 
of the local strong interaction dynamics. However, the 
presence of the electromagnetic field can affect the in- 
teractions among particles, which simultaneously carry 
electric and (possibly) color charges. 

Using relations ([5]) , the system ^ is reduced to a more 
familiar form: 



(6) 



(-Vr [/ + eE + evxB)Vp } f = Ccouifji, -fN) 

for particles of a charge e. 

The HSD transport model is based on Kadanoff-Baym 
equations for Green's function accounting for the first 
order gradient expansion of the Wigner transformed 
Kadanoff-Baym equation [37l - l39| . The hadronic mean 
field in Eq.© is C/ - ReiYT"^) /2po where E''^* denotes 
the retarded selfenergy. The change of the hadron mass 
in the HSD model results in an additional term ahead of 
Vp in (m (40j but it is ignored in our consideration which 
focuses on electromagnetic effects. 



Transport equations for a strongly interacting particle 
with electric charge e (EqH]) are supplemented by the 
electromagnetic field equations 



V x E = - 



1 9B 



V-B = 



(7) 



The general solution of the wave equations (O with 
the charge distribution p{v, t) = en and electric current 
j(r, t) = ev are 



for the electromagnetic scalar potential <I'(r,t) and 



A(r,i) 



47r 



i{r',t')5{t-t' -\v- 



(Pr'dt' (9) 



for the vector potential. For a moving point-like charge 
one gets 

p(r, t) = e5{v- r(t)), j(r, t) = e v(i) 5{v - r(t)) (10) 
and, after integration of Eq.® using 



/OO / 
j{x) 5{f{x))dx^Y.^^ 



we obtain: 



$(r,i) 



ineo ^ R{t'Mt',) 
with the definitions 

<t',) = 1 



(11) 



(12) 



(13) 



In ([T3| the times t'^ are solutions of the retardation equa- 
tion 



By analogy, for the vector potential we get: 



A(r,t) 



47r 



E 



(14) 



(15) 



Thus, Eqs.® and (jH]) lead to the retarded Licnard- 
Wiechcrt potentials ((T^ and ((T5|) acting at the point 
R = r — r' at the moment t. The electromagnetic po- 
tentials $(r,i) and A{r,t) are generated by every mov- 
ing charged particle and describe generally the elastic 
Coulomb scattering as well as inelastic bremsstrahlung 
processes. Below we set £o = no =1 in Eas. ((T2|) and 
respectively. The retarded electric and magnetic 
fields can be derived now from Eqs. ((T^ and ([T51) using 
([5]) as follows: 
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E(r,t) 



e / n — b/c — [(n • v)n — v]/i? 



47r V K2ci? 
e / (-v(t')/c + n(t'))(wVc - n • v - i?/c(n • b)) 

4^ I K3ci?2 



(16) 



ret 



B(r,t) 



e /vxn b(t') X n(i') +v(t') X [(n- v)n- v]/i?, 
4^ Ki?2 + ^ K^cR ■ 

e / (v(t') X n(i'))(wVc - n • v - i?/c(n • b)) 
47r \ K^ci?^ 



(17) 



with the acceleration b = d/dt'v and the unit vector 
n — R/R. After simphfication and neglecting the accel- 
eration b we arrive at 

sign{e) a R(t) (1 - v'^/c^) 



eE 
eB 



[{R{t) ■ v/c)2 + i?2(t)(l - w2/c2)]3/2 ' 

sign{e) a [v x R(i)] (1 — v^/c^) 

C [(R(i) • V/C)2 + i?2(t)(l - z;2/c2)]3/2 



(18) 
(19) 



where in the left-hand side an additional charge e is in- 
troduced to get the electromagnetic fine-structure con- 
stant a ~ e2/47r = 1/137 in the right-hand site of these 
equations. The last expression is reduced to the familiar 
form of the retarded Lienard-Wiechert equation for the 
magnetic field of a moving charge 



B(r,i) 



e 

47r 



X R] 



(l-t;2/c2) 



ci?3 [1 - (w/c)2sin2 0j^^]3/2' 



(20) 



with R ~ R(i), (j)Rv is the angle between R(i) and v. 

The set of transport equations in the following 
is solved in a quasiparticle approximation by using the 
Monte-Carlo parallel ensemble method. To find the elec- 
tromagnetic field, a space grid is used. The quasiparticle 
propagation in the electromagnetic field is calculated ac- 
cording to Eq.® as 



dp 



X B 



(21) 



The change of the electromagnetic energy is (e/c)(v • E) 
i.e. the magnetic field does not change the quasiparticle 
energy. To avoid singularities and self-interaction effects, 
particles within a given Lorentz-contracted cell are ex- 
cluded from the field calculation. 



III. RESULTS 

The scheme described above for the computation of the 
electromagnetic field is applied here to ultra-rclativistic 
heavy-ion collisions. However, for transparency, we shall 
start with the magnetic field created by a single freely 



moving charge. As seen in Fig. [U the charge creates 
a cylindrically symmetric field with the symmetry axis 
along the direction of motion. If one follows the mag- 
netic field direction, it appears to be torqued around the 
direction of motion. Therefore, the magnetic field on the 
left and right sides with respect to the moving charge has 
opposite signs, resulting in some maximum and minimum 
of the magnetic field at a given instant of time. The op- 
posite field signs directly follow from Eq. (|20|) if one takes 
into account that the vector R fEa. ([T^ ') in this situation 
has opposite signs. 
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FIG. 1: (Color online) Snapshot of the By distribution for the 
magnetic field and its projection on the {z — x) plane for a 
single charge moving along the z axis. 

In a nuclear collision, the magnetic field will be a 
superposition of solenoidal fields from different moving 
charges. The collision geometry for a peripheral collision 
is shown in Fig. [2] in the transverse plane. The over- 
lapping strongly interacting region (participants) has an 
"almond" -like shape. The nuclear region outside this al- 
mond (shaded in Fig. [2]) corresponds to spectator matter 
which is the dominant source of the electromagnetic field 
at the very beginning of the nuclear collision. Note that 
in the HSD code the particles are subdivided into target 
and projectile spectators and participants not geomet- 
rically but dynamically: spectators are nucleons which 
suffered yet no collision. 
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FIG. 2: (Color online) The transverse plane of a noncentral 
heavy-ion collision. The impact parameter of the collision is 
denoted by b. The magnetic field is plotted by the dashed 
lines. 



A. Space-time evolution of the magnetic field 



The time evolution of eBy{x, y — 0, z) for Au+Au col- 
lisions for the colliding energy ^/SNN =200 GeV at the 
impact parameter b =10 fm is shown in Fig. [31 If the 
impact parameter direction is taken as the x axis (as in 
the present calculations), then the magnetic field will be 
directed along the y-axis perpendicularly to the reaction 
plane {z — x). The geometry of the colliding system at 
the moment considered is demonstrated by points in the 
(z — x) plain where every point corresponds to a spectator 
nucleon. It is seen that the largest values of eBy ^ 5m^ 
are reached in the beginning of a collision for a very 



AuAu, Vs^jv = 200 GeV, b=10fm, t=0.01 fm/c AuAu, Vs^^=200GeV, b=10fm, t=0.05 fm/c 




AuAu, VS™=200GeV, b=10fm, t=0.2 fm/c AuAu, Vs»,», = 200 GeV, b=10fm, t=0.5 fm/c 




FIG. 3: (Color online) Time dependence of the spatial distribution of the magnetic field By at times t created in Au-|-Au 
(-^5=200 GeV) collisions with the impact parameter b —10 fm. The location of spectator protons is shown by dots in the 
{x — 2)-plane. The level By = and the projection of its location on the {x — z) plane are shown by the solid lines. 



short time corresponding to the maximal overlap of the 
colliding ions. Note that this is an extremely high mag- 
netic field, since « lO^^gauss. The first panel in 
Fig. [3] is taken at a very early compression stage with 



t =0.01 fm/c. The time t =0.05 fm/c is close to the 
maximal overlapping and the magnetic field here is max- 
imal. Then, the system expands (note the different z- 
scales in different panels of Fig. [3]) and the magnetic field 
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decreases. For b =0 the overlapping time is maximal and 
roughly given by 2R/^c which for our case is about 0.15 
fm/c. For peripheral collisions this time is even shorter. 

Globally, the spatial distribution of the magnetic field 
is evidently inhomogeneous and Lorentz-contracted along 
the z-axis. At the compression stage there is a single 
maximum which in the expansion stage is splitted into 
two parts associated with the spectators. In the trans- 
verse direction the bulk magnetic field is limited by two 
minima coming from the torqued structure of the single- 
charge field (see Fig. [T]). 

The possibility of attaining extremely high magnetic 
fields in heavy-ion collisions was pointed out 30 years 
ago [m but there have been only two real attempts to 
estimate the magnetic field for relativistic heavy-ion colli- 
sions [111, [2^ . In Ref. [2l[ the colliding ions were treated 
as infinitely thin layers (pancake- like) , and the results 
in the center of a Au-Au collision eBy{0,0, z) could be 
presented in a semi-analytical form. In Fig. 2] these esti- 
mates are confronted with our results. It is clearly seen 
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FIG. 4: (Color online) Time dependence of the |ei3| field in 
the center of the nuclear overlap region for Au+Au(y^=200 
GeV) collisions from the HSD calculations. The dotted and 
dot-dashed curves are from Ref. [SJ at the impact parameters 
b =8 and 12 fm, respectively. 

that the magnetic field in our transport model for 6 =10 
fm is lower than the estimate from Ref. [2l| for both 
6 =12 and 8 fm. This difference originates mainly from 
the fact that to simulate rapidity degradation of pancake- 
like nuclei, a heuristic function was assumed with making 
no difference between surviving baryons and new created 
particles [2lj whereas in our case the dynamical hadron- 
string model is used for both primary and subsequent 
interactions while keeping electric and baryonic charges 
and energy- momentum conservation [36j . The approxi- 
mation of Ref. [Ill is reasonable for first collisions but 
gets progressively worse with interaction time as seen in 
Fig. m The difference in the shape of the time depen- 
dence of the magnetic field for early times is due to ne- 
glecting the finite size of the colliding nuclei in Ref. [2l| . 
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FIG. 5: (Color online) The magnetic field evolution at the 
point X for y —0. 



AuAu, = 200 GeV, b=10 fm 
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FIG. 6: (Color online) Time evolution of the magnetic field 
at the point y for the central overlap point x =0. 



Also, in our treatment the self-interaction is excluded 
for charges within the Lorentz-contracted hadron volume. 
Our consideration treats more accurately the retardation 
effect discussed above which constrains the contributions 
to the given point from some charges. It is especially 
important for the field contribution from participants. 

It is of interest to note that in our transport model, 
the spectator contribution to the magnetic field is prac- 
tically vanishing at i wl fm/c (see Fig. [3]). In subsequent 
times the magnetic field eBy is formed essentially due 
to produced participants with roughly equal number of 
negative and positive charges which approximately com- 
pensate each other. The visible effect in our approach 
is by an order of magnitude lower than that in the es- 
timate [2l| which demonstrates the essential role of the 
retardation in this interaction phase. 

Furthermore, the magnetic field distribution in (26j 
is calculated within the UrQMD model and the back 
reaction of the field on particle propagation is disre- 
garded. Nevertheless, our results are quite close to those 
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of Rcf. [1^. 

In Fig. [51 the magnetic field evolution eBy{x, y = 0,z) 
is shown as a function of the transverse coordinate x. 
Practically the difference between results for x =0, 1, 2 
fm is less than 20% except for the boundary of the overlap 
region corresponding to a; « b/2 ~5 fm. One thus may 
conclude that the magnetic field is rather homogeneous 
in the transverse direction. 

The magnetic field component By{x = 0, y, z) along 
the largest axis y of the "almond" (see Fig.[2|) is presented 
in Fig. ini for different times. The similarity of all curves 



for y <4 demonstrates a high homogeneity of the created 
field By. It is of interest that this field stays almost 
constant during At ^0.1 fm/c. 



B. Energy density and its correlation with By 

Along with a high magnetic field, the presence of a 
quark-gluon phase is a necessary condition for a manifes- 
tation of the chiral magnetic effect according to 



AuAu, V5»,„=200GeV, b=10 fm, t=0.01 fm/c 



0.3 
0.2 
_ 0.1 

t§ 

N 

-0.1 
-0.2 
-0.3 



-1 1 1 1 r- 



_l I I I L_ 



10 „ 

> 

CD 

1 a 



0.1 



-15 -10 -5 5 
X [fm] 



10 15 



AuAu, V5;^,;v=200GeV, b=10 fm, t=0.2 fm/c 



0.3 
0.2 
0.1 


-0.1 
-0.2 
-0.3 



— I ' ' 



44 



1 



_l I I L_ 



10 



0.1 



AuAu, V5»,„=200GeV, b=10 fm, t=0.05 fm/c 



0.3 
0.2 
0.1 
I 

N 

-0.1 
-0.2 
-0.3 



1 1 1 1 1 

. .eBy/m^^ 


1 1 

,-' 




r3 - 






- '-''o 

1 1 1 1 1 





10 _ 

> 

CD 

1 o. 



L-l 0.1 



-15 -10 -5 5 
X [fm] 



10 15 



AuAu, V5^;v=200GeV, b=10 fm, t=0.5 fm/c 







0.6 






0.4 


"'p 




0.2 


eV/fi 


[fm] 





O 


N 


-0.2 






-0.4 
-0.6 



1 1 1 


1 1 1 1 










e By/m^^ J 


I 


'.0 


1 

1 1 1 


-1, , - 



10 ^ 

> 

<D 

1 a 



0.1 



-15 -10 -5 5 
X [fm] 



10 15 



-15 -10 -5 5 
X [fm] 



10 15 



FIG. 7: (Color online) Spatial correlations in location of the magnetic field and energy density of participants produced in 
Au+Au(y^ =200 GeV) collisions witli the impact parameter 6 =10 fm at the times t. The levels of the magnetic field are 
plotted by dashed lines whereas areas with different energy densities are displayed in color. 



Refs. [isl [l9l - l23i . The phase structure of excited matter 
is essentially defined by the energy density (cf. Ref. ji^). 
One can expect that for energy densities e >1 Gev/fm^ 
the system is in a deconfined phase. The evolution of the 
energy density of created particles is presented in Fig. [T] 
Here the maximal energy density (in the center of the 
colliding system) is e > 50 GeV/fm^ at the moment 
of maximal overlap of the nuclei. When the system ex- 
pands, it takes a sausage-like shape (or dumb-bell shaped 
if the energy density values are taken into consideration 
additionally) and the energy density drops fast. But even 
at the time t ~ 0.5 fm/c (last panel in Fig. [7]) the local 



energy density is seen to be above an effective threshold 
of a quark-gluon phase transition e >1 GeN / fm? . Dif- 
ferent levels of the magnetic field strength are plotted in 
the same figure. It is clearly seen that the location of 
the maximum energy density correlates with that for the 
magnetic field. 

The variation of the energy density distribution with 
the transverse coordinate x is shown in Fig. [51 Here 
the plotted values of e correspond to averages within the 
Lorentz-contracted cylinder with \z\ < 5/7 fm and ra- 
dius R =1 fm centered at the point x. One can see that 
the energy density changes more strongly in x than the 
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AuAu, V5'„«,= 200GeV, b=10fm 




FIG. 8: (Color online) The average energy density in the 
Lorentz-contracted cylinder of the radius R = 1 fm and 
\z\ <5/7 fm with the z-axis passing through the point x. 



transverse momentum < ptr > during the whole evolu- 
tion even if the fastest pions escape the finite cylinder 
volume by the time t ^ 0.2 fm/c. The ratio < pi > / < 
Ptr > is very large and does not correspond to that in 
equilibrium. The sharp decrease of < p; > in time is due 
to fast mesons streaming-out from the finite volume due 
to the longitudinal expansion of the system rather than 
to equilibration. Thus, in the early stage t <0.2 fm/c 
the Au-f-Au system with high energy density is far from 
equilibrium. On the one hand, this fact is not astonishing 
since this stage should be treated in terms of quarks and 
gluons rather than on a hadronic level (cf. correspond- 
ing PHSD studies On the other hand, there is a 
general consensus that local equilibrium hydrodynamics 
can be applied to heavy-ion reactions at energies cur- 
rently available the BNL Relativistic Heavy Ion Collider 
(RHIC) only for t > to «0.5 fm/c [H,!!!. 
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FIG. 9: (Color online) The average longitudinal and trans- 
verse momenta of all mesons and pions in the Lorentz- 
contracted cylinder described above. 



magnetic field (note the logarithmic scale in Fig. [5]) . In 
particular, the maximal e decreases by a factor ~20 when 
one proceeds from a; = to a; =3 fm and close to the 
spectator-participant boundary (at x ~3 fm) the energy 
density very quickly (within roughly ^0.3 fm/c) drops 
below the effective threshold for deconfinement e ^1 
Gev/fm'^. 

One should note that the energy density should be cal- 
culated in the rest system. The choice of a symmetric 
position of the cylinder volume with respect to the z =0 
plane essentially leads to an approximately vanishing to- 
tal momentum of particles inside this volume. The time 
averaged 7-factor of particles in this cylinder in the cm. 
system is < 7 >~T1. Note, however, that the created 
particles are not in local equilibrium! 

In Fig. [S] the evolution of the average longitudinal 
< pi > and transverse < ptr > momentum is shown for 
mesons in the same cylinder. All mesons keep a constant 



C. Electric field evolution 

The background electric field, being orthogonal to the 
magnetic one, is directed along the x axis. The evo- 
lution of the eEx field for peripheral (6 =10 fm) colli- 
sions of Au-|-Au at the top RHIC energy is presented in 
Fig. [TOl Similar to the case of the magnetic field, the 
eEx{x,y = 0,z) distribution is also inhomogeneous and 
closely correlates with geometry while the field strength 
looks "hedgehog" shaped. When the two nuclei collide, 
the electric fields in the overlap region significantly com- 
pensate each other, and the electric field E in the target 
and projectile spectator parts have opposite signs. As a 
result, the locations of the maximum/minimum are not 
in the central point of the overlap region - as they are 
for the magnetic field - but shifted slightly outside. The 
maximum of the electric field can be quite large. All these 
features are seen explicitly in Fig. [TT] where the tempo- 
ral evolutions of eEx{x, 0, 0) and eEy{x, 0, 0) are given for 
different values of the transverse coordinate x. Due to de- 
structive interference or the "hedgehog" effect, the elec- 
tric field in the central part of the overlap region {x 
fm) is consistent with zero apart from a short period just 
before reaching maximal overlap. For x ~l-3 fm the elec- 
tric field has a distinct maximum of eEx/m^ w(0.5-0.6) 
which is only by a factor of about 10 less than the max- 
imal magnetic field eBy (cf. Fig. [S]). But when moving 
farther away from the center of the overlap region the 
eEx component drops sharply and then becomes neg- 
ative. The eEy field component is quite small for the 
central part of the overlap region and increases slightly 
for larger x. Note that the electric field is negligible for 
t >0.15 fm/c. 
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FIG. 10: (Color online) Evolution of the x- and y-components of the electric field at incoming and maximal overlap in 
Au+Au(-ys= 200 GeV) collisions at the impact parameter h =10 fm. The eEx =const levels and spectator points are shown 
in the projection on the (x — z) plane. 
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FIG. 11: (Color online) Evolution of the x (left panel)- and y (right panel)-components of the electric field for different values 
of the transverse coordinate x. 



D. Scalar product of (E • B) ternal magnetic field. A magnetic field B can also com- 

bine with an electric field E to form the Lorentz invariant 
(E • B) which changes the sign under a parity transfor- 
Since the magnetic field is odd under time reversal (or mation. In the normal QCD vacuum with its sponta- 
equivalently, under the combined charge conjugation and ncously broken chiral symmetry the leading interaction 
parity CV transformation), the time reversal symmetry involves the invariant (E • B) which enters e.g. into the 
of a quantum system is broken in the presence of an ex- 
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FIG. 12: (Color online) Space-time evolution of of the scalar product of electric and magnetic field (E-B) for a Au+Au reaction 
with the impact parameter b =10 fm at -^5=200 GeV. 



matrix element that mediates the two-photon decay ol 
the neutral pseudoscalar mesons. In the deconiined chi- 
rally symmetric phase of QCD, the leading interaction 
term is proportional to aas^E ■ B)(E'* • B"^), where E"^ 
and B"^ denote the chromoelectric and chromomagnetic 
fields, respectively, and a and as are the electromagnetic 
and strong QCD couplings. Both interactions are closely 
related to the electromagnetic axial anomaly, which in 
turn relates the divergence of the isovector axial current 
to the pseudoscalar invariant of the electromagnetic field 
(see Ref. [iH). The evolution of the electromagnetic in- 
variant - produced in Au+Au (•y/i=200 GeV) collisions 
at the impact parameter b =10 fm - is shown in Fig. 1121 



As seen the electromagnetic invariant (E • B) is non- 
zero only in the initial time t <0.5 fm/c where the (E-B) 
distribution is quite irregular and its non-zero values cor- 
relate well with the location of the overlap region. For 
longer times this electromagnetic invariant vanishes as 
follows from the electric field space-time distributions 
(cf. Figs, [in] and [TT|) . Note that the quantities plot- 
ted in Fig. [T2] are dimensionless and the scaling fac- 
tor [Gey^] is quite small. In a topological domain 
the chromoelectric fields (E° • B°) y^O. But here non- 
vanishing values of the electromagnetic invariant are due 
to highly inhomogeneous distributions of electromagnetic 
fields E and B. 



E. Impact parameter dependence 

As noted above, the electromagnetic field is formed 
predominantly by spectators during the passage time of 
the two colliding nuclei. Since the number of specta- 
tors increases with the impact parameter b, the magnetic 
field should also increase for more peripheral collisions. 
Indeed, as seen in Fig. [131 the magnetic field decreases 
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FIG. 13: (Color online) Impact parameter dependence of the 
magnetic field in Au-|-Au collisions at ^/snn =200 GeV. 
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FIG. 14: (Color online) Impact parameter dependence of the 
eEx{x = 0,0,2) (top panel)- and eEy{x = 0, 0, z) (bottom 
panel)- components of the electric field in Au-|-Au collisions 
at ^/i^ =200 GeV. 



above the e£'j,-component; both are constrained to the 
time interval < t < 0.2 fm/c (Fig. [T4|) . Irregularities 
in these distributions are due to the "hedgehog" effect 
mentioned above. 

In contrast, the impact dependence of the energy den- 
sity £ should correlate with the number of participants 
and therefore reach a maximum in central collisions. The 
^-dependence of the e temporal distributions is quite 
weak (see Fig. [15]), and it decreases (by a factor of ^2) 
only for far peripheral collisions, 6 ~10 fm. 

Note that in the {x— z) plane, the location of the max- 
imal values of the magnetic field and the energy density 
correlate with each other (see Fig. [7]). 



F. Dependence on collision energy 

In principle, the collision energy dependence of elec- 
tric and magnetic fields is given by Eqs. (|16p and (|17p . 
respectively, since the fields predominantly emerge from 
spectators moving with the initial velocity v. These ex- 
pressions involve velocity-dependent factors v and 1 — u^, 
which vanish in the limiting cases v —^0 or w — >1. In any 
of these limiting cases, the denominators stay finite and 
the energy dependence is given by the numerator. So, the 
electric field approaches zero in the ultrarelativistic limit 
V — >-l; the magnetic field vanishes in both limits v — >-0 
and V ^1. As follows from Fig. [1^1 the maximal strength 



gradually with decreasing h. When changing the impact 
parameter from b =10 to 2 fm, the maximal eBy{x, 0, 0) 
decreases by a factor of 5. For t > 0.3 fm/c the magnetic 
field goes down almost exponentially (cf. with Fig. |4]). 
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FIG. 15: (Color online) Impact parameter dependence of the 
average energy density in the Lorentz-contracted cylinder of 
radius i? = 1 fm and \z\ < 5/7 fm with the z-axis passing 
through the point 2; = in Au+Au collisions at ^/sNN =200 
GeV. 

For the electric field, the ei?^-componcnt is only slightly 
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FIG. 16: (Color online) Collision energy dependence of the 
magnetic field evolution in the center point of the overlap 
region. 

of the magnetic field ei?j,(0, 0, 0) decreases roughly pro- 
portionally to ^sjvAf and at the top energy available at 
the CERN Super Proton Synchrotron (SPS) (y/s^N «18 
GeV), it is only about ~ 0.4/m^ which appears to be too 
low to search for the Chiral Magnetic Effect 27|. One 
should keep in mind that the chiral magnetic effect de- 
pends not only on the magnitude of the magnetic field but 
also on the time of the system within the high magnetic 
field. As seen in Fig. [121 with decreasing collision energy 
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FIG. 17: (Color online) Collision energy dependence of the x 
and y component of the electric field. 
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FIG. 18: (Color online) Collision energy dependence of the 
average energy density in the Lorentz-contracted cylinder of 
radius 7? = 1 fm and \z\ < 5/y fm with the z-axis passing 
through the point a:: = 0. 



the width of the time distribution of By becomes wider 
because it is defined essentially by the Lorentz-contracted 
overlap region ^ Ar/7 — / y/s^N. The RHIC 

energy is not high enough to see the limiting case corre- 
sponding to (1 — w^) — )-0. 



According to our expectation, the maximal strength 
of the eEx{Q,0,0) and eE'j, (0,0,0) components decrease 
for lower collision energy as demonstrated in Fig. [T71 It 
is clearly seen by comparing the results for ^/snn =200 
and 18 GeV. Note that the eEy component is directed 
opposite to the direction of the magnetic field eBy and 
will act against "the electric charge separation effect" j^^, 

MM- 

The maximal energy density drops down by two orders 
of magnitude when going down from the RHIC to the 
SPS energy (see Fig. [75]). During the time interval 0.5< 
t <2.5 the energy density e{t) evolution for all energies 
practically coincide changing from e '-^lO to '-^l GcV/fm^. 



IV. OBSERVABLES AND ELECTRIC CHARGE 
SEPARATION 

The HSD model quite successfully describes many ob- 
servables in a large range of the collision energy [la 1131 ■ 
Here we investigate to what extent the electromagnetic 
field - incorporated in the HSD approach - will affect 
some observables. We shall limit ourselves to Au-|-Au col- 
lisions at -y/J/vTV =200 GeV and impact parameter b =10 
fm. In this case we arc not able to restrict our calcula- 
tions to times t <3 fm/c as above but have to calculate 
the whole nuclear interaction including the decays of res- 
onances at least up to times of 50 fm/c. 

The HSD results for the versions without and with elec- 
tromagnetic field - taking into account the back reaction 
of the clcctrodynamic field on the particle propagation - 
are presented in Fig. (THl With a high degree of accuracy, 
we see no difference between these two versions in the 
transverse mass rrit and rapidity y. 

In Fig.[20]thc transverse momentum dependence of the 
elliptic fiow of charged pions is compared for two versions 
(with and without field) of the HSD model. We do not 
observe any significant difference between the two cases. 
Slight differences are seen in the range oi pt ^1 GeV/c 
but certainly it can not be considered as significant. Note 
that generally the HSD model underestimates the ellip- 
tic flow, but an inclusion of partonic degrees of freedom 
within the PHSD approach allows it to describe perf ectly 
well the Pt dependence of V2 at the top RHIC energy [i^l ■ 

As a signal of possible CV violations in relativistic 
heavy-ion collisions, it was proposed in Ref. (46j to mea- 
sure the two-particle angular correlation 



(cos(0„ +(f>(i- 2*flp)) 



(22) 



where 5'/fp is the azimuthal angle of the reaction plane 
defined by the beam axis and the line joining the cen- 
ters of colhding nuclei (see FigH]). The correlator ([^ is 
calculated on the event- by- event basis with subsequent 
averaging over the whole event ensemble. The experi- 
mental data from the STAR Collaboration [2^, [2^ and 
the results of HSD calculations are presented in Fig. [5T] 
The experimental acceptance \r]\ < 1 and 0.15 < pt < 2 
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FIG. 19: (Color online) Transverse mass and rapidity distri- 
butions of charged pions produced in Au+Au {^/snn =200 
GeV) collisions at b =10 fm. The results calculated with and 
without electromagnetic field are plotted by the dotted and 
solid lines, respectively. 



tions. Note that the theoretical reaction plane is fixed 
exactly by the initial conditions and therefore is not de- 
fined by a correlation with a third charged particle as in 
the experiment [1^ [2^ . The error bars plotted in Figl^T] 
show the statistical errors. The number of events evalu- 
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FIG. 20: (Golor online) The transverse momentum depen- 
dence of the elliptic flow for Au+Au {^/snn =200 GeV) col- 
lisions at b =10 fm. 



GeV have been also incorporated in theoretical calcula- 



FIG. 21: (Color online) Azimuthal correlation in the trans- 
verse plane as a function of centrality for like and unlike 
charged pions from Au+Au {^snn =200 GeV) collisions. 
The experimental points - connected by lines - are taken 
from m,^- 



ated with (without) field for the most crucial centralities 
0.7 and 0.55 are e.S-lO^ and 2.2-104 (8.4-104 and 5.4-10^) 
for same charge pion pairs. The computational time for 
one event with accounting for the electromagnetic field 
is by a factor of about 30 longer than that for the case 
without field. 

The expected CME stems from the interplay of topo- 
logical effects of the excited vacuum and the chiral 
anomaly in the presence of a strong magnetic field (isl . 
[20l - [23l |. One can see that the calculated background - 
taking into account hadron string interaction dynamics 
and evolution of the electromagnetic field - is not able to 
describe the measured distribution especially for pions of 
the same charge. One should mention that our results 
arc rather close to the background estimates within the 
UrQMD model in the experimental works [2^, [1^ . 
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FIG. 22: (Color online) Projections of azimuthal correlations on the in- and out-reaction plane as a function of centrality for 
like and unlike charge pious from Au+Au (y^ijvjv =200 GeV) collisions. The experimental points - connected by lines - are 
from Ref. S. 



The two particle correlation (p2|) can be decomposed 
in "in-plane" and "out-of-plane" components [52j | 



(C0s((?!)a + 0^)) = (cOs((^a)cOs(0/3)) 

- (sin((/)Q)sin((?!)^)) . 



(23) 



Following Ref. [34] in Fig.[22]thcse components are pre- 
sented for the same (+,+),(—,—) and opposite (+,— ) 
charged pion pairs. First, there is no difference for HSD 
results without field (top panels in Fig. [22])and with the 
electromagnetic field (bottom panels). To be more spe- 
cific, we will discuss below the results with the electro- 
magnetic field included. Second, the calculated differ- 
ence between the in-plane (cosme-term) and out-of-plane 
(sine-term) components is small in the same and oppo- 
site charge cases. 

Since the observed correlation (|23p is the difference of 
these two terms, the calculated correlation is small as 
well. Furthermore, for the same charge pairs the mea- 
sured sine term is essentially zero while the cosine term 
is finite. This implies that the observed correlations are 
in-plane rather than out-of-plane, as expected. It is of 
interest that the measured and calculated cosine terms 
coincide with each other for centralities <0.55. As was 



noted in Ref. |32| . the zero sine component is contrary 



to the expectation from the CME, which for the same 
charge correlation results in an out-of-plane correlation. 
In the HSD model the sine term is not zero but negative. 
This is not a surprise because the induced chromoelectric 
field parallel to the out-of-plane By is not included into 
our calculations, but there is a nonzero electric field com- 
ponent Ey (see above). Furthermore, we see that for op- 
posite charge pairs the sine and cosine correlation terms 
are virtually identical, which, according to Refs. [32l.[33|. 
is hard to reconcile with a sizable elliptic flow in these 
collisions. However, the centrality distributions of op- 
posite charge pions exhibit contrary trends: the STAR 
measurement is positive and decreases but the HSD re- 
sult is negative and increases toward central collisions 
where all components of angular correlations wO. It is 
noteworthy that the UrQMD model shows quite close re- 
sults. Indeed, (cos(0q, — (/)^)) is just the sum of cosine and 
sine terms. So, summing the two opposite charge curves 
in Fig. [21] we reproduce the UrQMD results presented in 
Fig. 5 of Ref. Hi. 

Recently there were proposals to explain the observed 
CME effect without reference to the local parity viola- 
tion. Different background mechanisms of the azimuthal 
correlation are considered: cluster decays, local trans- 
verse momentum conservation and local electric charge 
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conservation [30l - l33| . These mechanisms may contribute 
to the effect under discussion, but their simphfied esti- 
mates made in these papers are not able to describe the 
STAR measurements. GeneraUy speaking, aU these ef- 
fects, such as the decay of resonances including heavy 
ones and the exact conservation of electric charge and 
energy-momentum, arc involved in our transport model 
but they do not help. However, there is a conflicting 
point here. The HSD model treats the system evolution 
for Au-|-Au(200 GeV) in terms of hadrons and strings 
while the decisive phenomena occur for times t <0.3 
fm/c. This is definitely in a non-equilibrium quark- 
giuon state. On the other hand, as was shown in the 
multi-phase model [i^ , the charge separation can be sig- 
nificantly reduced by the evolution of the quark-gluon 
plasma produced in relativistic heavy-ion collisions and 
by the subsequent hadronization process. 

The conservation of the transverse momentum is a 
possible source of azimutlial correlations which was sug- 
gested to be a significant contribution to the measured 
observable (Eq.^H) dH-dl]). Using Eq.l^ and the con- 
servation of the transverse momenta, it was shown that 
roughly PIH) 

(cos(</)„ + 00-21'^p)) ~ , (24) 

where V2 is the elliptic fiow coefficient measured for all 
produced particles and N is the total number of all pro- 
duced particles. It is definitely a qualitative result, but 
it demonstrates a close relation of the observed charge 
separation effect with the elliptic fiow. Equation ([24|) 
was specified more accurately in Ref. [11]) and the corre- 
lator ([22)) was estimated under reasonable assumptions. 
It was concluded that transverse momentum conserva- 
tion alone is not sufficient for explaining the STAR data. 
One should add that this issue has been considered in the 
multi-phase model [1^ as well, with the conclusion that 
the charge particle separation leads to a modification of 
the relation between the charge azimuthal correlation and 
the elliptic fiow that is expected from transverse momen- 
tum conservation (Eq. [24)) . An essential point from this 
discussion is that the charge separation effect is roughly 
proportional to the elliptic flow; however, this quantity 
is underestimated in the HSD model resulting in a small 
value of the correlator. We hope that future calculations 
within the PHSD model [i^ might provide azimuthal 
correlations that will be closer to the measured data. At 
present the questions noted above regarding the experi- 
mental observations have no simple explanation. 

V. SUMMARY AND OUTLOOK 

We have extended the hadron string dynamics model 
to describe the formation of the retarded electromagnetic 
field, its evolution during a nuclear collision and the ef- 
fect of this field on particle propagation. The case of 
the Au-|-Au collision at ^/snn for b =10 fm is consid- 
ered in great detail. It is quite important to understand 



the interplay of strong and electromagnetic interactions 
in this case since it provides a point which is decisive in 
the CME measurements at RHIC as a function of cen- 
trality [2^, [2^ . It is shown that the most intensive mag- 
netic field oriented perpendicularly to the reaction plane 
is formed during the time when the Lorentz-contracted 
nuclei arc passing through each other, t <0.2 fm/c. The 
maximal strength of the magnetic field here attains very 
high values, eBy/m^ ^5. This magnetic field strength 
is higher by about 3-4 orders of magnitude than that 
at the surface of the magnetar [49j which in turn is only 
slightly above the field strength in the star core [s^l . Still 
larger magnetic fields up to i? ^ lO^'^ G might have ap- 
peared in the early universe [Elj . It is impossible to make 
steady fields stronger than 4.5 • 10^ gauss in the labora- 
tory because the magnetic stresses of such fields exceed 
the tensile strength of terrestrial materials. 

This maximal strength of the magnetic field is created 
predominantly by spectators. When target and projec- 
tile remainders are separated, the spectator contribution 
goes sharply down and for t ^1 fm/c decreases by more 
than three orders of magnitude. In subsequent times the 
participants come into the game, but their contributions 
are small due to the mutual compensation of approxi- 
mately equal number of positive and negative charges as 
well as to the suppressive role of the relativistic retarda- 
tion effect. 

The general pattern of the magnetic field is highly in- 
homogeneous. However, in the "almond" transverse area 
(besides the boundary region) the field distributions in 
z or time look very similar, which allows us to use some 
simplifying assumptions in phenomenological CME esti- 
mates i26|,i27|. 

The important accompanying quantity is the energy 
density e of the created particles. Its space-time dis- 
tributions have been presented. It was shown that the 
location of maxima in the field strength eBy and the en- 
ergy density e nicely correlate with each other. Thus, it 
is a necessary condition for a realization of the CME. 

In the early time moments the created electric field is 
perpendicular to the magnetic one and has a dominant x- 
component. In contrast to the eBy distribution, the eE^ 
distribution in the (x — z) space plane has a minimum in 
the center of the overlap region due to the "hedgehog" 
field structure of an isolated electric charge. The maxi- 
mal strength of the electric field is by a factor about 5 
lower than that of the magnetic field. For t >0.20 fm/c 
the electric field is small and can be even neglected. 

The electromagnetic field is only moderately (within 
factor ^5) changed with impact parameter and colli- 
sion energy (for ^Js^n ~200-60 GcV) and strongly sup- 
pressed for t >0.2 fm/c. The scale of the energy den- 
sity change is much larger but up to times of a few fm/c 
£ >1 GeV/fm^ and thereby does not prevent quark-gluon 
plasma formation. Certainly, the issue of thermalization 
remains open in this consideration. The very small infiu- 
ence of the electromagnetic field on observables is due to 
rather large masses of quasiparticles and the point that 
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the system spends a very short time in a state with an 
extremely high electromagnetic field. 

The comparison of global observables calculated in the 
HSD model with and without an electromagnetic field re- 
veals no difference apart from the transverse momentum 
dependence of the elliptic flow where the model results 
slightly differ in the range pt ^1 GcV/c. 

Our analysis of the angular correlators - specific for 
the CME shows that the calculated HSD background is 
very small and not able to describe the STAR measure- 
ments [2^ [2^ . The consideration of in-plane and out-of- 
plane projection components of this correlator does not 
allow us to clarify the picture and rises new questions 
related to the experiment [13, • In this respect it is of 
great interest to include quark-gluon degrees of freedom 
directly in our approach. In particular, the partonic gen- 
eralization of the HSD model (PHSD) is highly suited 
for this aim. Another way to approach the CME is to sim- 



ulate an induced chromoclcctric field which is assumed to 
be the source of the observed electric charge separation 
of pion pairs relative to the reaction plane. Alternative 
mechanisms of charge azimuthal asymmetry should be 
also carefully studied. Certainly, measurements of the 
CME at other bombarding energies as well as a search 
for new observables are very important. 
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